Phasor analysis of atom diffraction from a rotated material grating 



m 
o 
o 

(N 
O 

Q 



^ ■ 

Oh: 



O ■ 

>v 
43 ■ 



> ■ 
■ 

oo ' 

o : 

(N ■ 

m : 
o . 

o : 

>>: 
43 : 

Oh. 



X 



Alexander D. Cronin, John D. Perreault 

University of Arizona, Tucson, Arizona 05721 
(Dated: February 2, 2008) 

The strength of an atom-surface interaction is determined by studying atom diffraction from a 
rotated material grating. A phasor diagram is developed to interpret why diffraction orders are never 
completely suppressed when a complex transmission function due to the van der Waals interaction is 
present. We also show that atom-surface interactions can produce asymmetric diffraction patterns. 
Our conceptual discussion is supported by experimental observations with a sodium atom beam. 
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Atom diffraction from a material grating has recently 
been used to measure the strength of van der Waals inter- 
actions between atoms and the grating. This was possible 
in 0, 0, 13 because atom-surface interactions modify the 
intensity in each diffraction order as a function of atomic 
velocity. In this paper we measure van der Waals inter- 
actions by studying diffraction as a function of incidence 
angle. We also demonstrate an asymmetric atom diffrac- 
tion pattern from a fabricated material grating. 

Historically, a graphical analysis in the complex plane 
has been useful to understand optical diffraction. This 
is especially true for the Fresnel integrals for which no 
closed form analytical solutions have been found, yet the 
Cornu spiral permits a physical interpretation [4j. We 
have adapted this approach to our current problem in 
atom optics. Even the far-held limit, van der Waals in- 
teractions modify atom diffraction such that no closed 
form analytical solution has been found. Hence we de- 
veloped a phasor diagram similar to the Cornu spiral to 
interpret our atom diffraction data. 

In particular, we prove that no combination of grat- 
ing geometry and van der Waals interaction strength can 
cause diffraction orders to disappear. This is proved us- 
ing a phasor diagram, and confirmed experimentally by 
rotating a diffraction grating about an axis parallel to the 
grating bars while measuring the flux in each diffraction 
order. 

In the standard optical theory of diffraction from a 
Ronchi- ruling, i.e. a square-wave absorbing grating, a 
missing order is obtained when 
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where n is the diffraction order, m is an integer greater 
than zero, and 77 is the open fraction defined as window 
size divided by grating period. For example, a 50% open 
fraction suppresses all the even orders. 

The origin of missing orders in the optical theory can 
be explained with a phasor plot often referred to as a 
vibration curve || and shown in Figure^,. In this ap- 
proach a complex amplitude-phase diagram is used to 
visualize the amplitude and phase of the held in each 
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FIG. 1: (a) Phasor diagrams for diffraction into orders 
n = 0, 1, 2, 3 from a grating with open fraction r\ — 0.48. Re- 
sultant vectors are drawn with dashed lines from tip to tail of 
each curve, (b) Van der Waals interactions modify the pha- 
sor diagram, and considerably increase the magnitude of the 
second order. 



diffraction order. Specifically, we consider the field fy n 
associated with the n th diffraction order in the Fraun- 
hauffer approximation due to a periodic array of aper- 
tures illuminated by a plane wave ^i nc - The slit width w 
and grating period d determine the diffraction envelope 
so the amplitude and phase of each order is given by: 
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where $ mc is the incident wave amplitude, £ is the spatial 
coordinate inside an aperture, and 4> n {C) = nn^ with 
k = being the grating wavenumber. 

This integral can be visualized by adding phasors of 
length d£ and phase </>„(£) in the complex plane. Curves 
generated this way are the real vs imaginary parts of 
the cumulative integral of ^> n given by Equation The 
magnitude and phase of a resultant vector (from start to 
end of such a curve) correspond to the complex amplitude 
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^ n as shown in Figure ^ A resultant of zero magnitude 
represents a missing order. 

Before including phase shifts due to atom-surface in- 
teractions, this integral can be computed analytically: 



= ^inc V sinc(n77) 
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For comparison with experiment, intensity is given by 
In = l*n| 2 - 

In the WKB approximation, to leading order in V(£), 
van der Waals interactions cause a phase shift 4> v dw(Q 
given by 
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where V(^) is the van der Waals potential for atoms in- 
teracting with a surface, I is the distance the atom prop- 
agates in the potential, v is the atomic velocity, and h 
is Planck's constant divided by 2ir. Near a surface the 
potential is known to be 



V(r) = 
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where r is the atom-surface distance For simplicity, 
as in 0,12,0] V{ r ) is approximated by equationEJinside a 
grating window with finite thickness I, and approximated 
by zero outside the grating window. 

Van der Waals interactions with the bars on either side 
of the slot thus modify the phase of each point on the 
phasor diagram. This phase is described by: 



4>n{£) = <f>n(€) + <t>vdw(£) + ^offset 
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where 4> n (£) is due to the diffraction angle, <fi v dw(Q is 
due to van der Waals interactions with surfaces located 
at £ = ±|r, and (fro ff set is a constant added for reasons 
that will be discussed: 
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Examples of phasor diagrams modified by atom-surface 
interactions are shown in Figure and for a range of 
van der Waals interaction strengths are shown again in 
Figure □ 

The constant <fioffset is not physical but is chosen to 
be 4>offset = — 4>vdw(Si) as a convenience to compare the 
shape of the phasor curves for arbitrary van der Waals 
constant C3 . This phase offset does not change the norm 
of the resultant amplitudes, hence it will not affect the 
proof regarding missing orders. It simply rotates the spi- 
ral such that </>„(0) = independent of C3. 

For C3 = the phasor curves are simple circles because 
the curvature = nK is constant. The phase <j) n {C) 
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FIG. 2: Phasor diagrams for diffraction from a grating 
with open fraction r\ = 0.48 shown for diffraction order 
n = 2 given different values of the van der Waals coefficient 
C 3 = 0, 1, 10, 100 eVA 3 . The resultant vectors from tip to tail 
are drawn with dashed lines. 



spans an angle unw, and the arc length of the curve is 
given by the window size w. Thus the radius of the circle 
is p = (nn)^ 1 and it is centered at the location ip. When 
the curve is an integral number of full circles, the end- 
points overlap and the resultant field has zero magnitude. 
This corresponds to a missing order. 

The additional phase due to van der Waals interactions 
makes the phasor curve deviate from a circle. One end 
of the spiral will always be inside the circle defined by p 
and the other end must be outside. This is true because 
the curvature -^{4> n + 4'vdw) increases monotonically as £ 
goes from to w/2 and decreases monotonically as £ goes 
from to -w/2 (and the curvature equals nn at £ = 0). 
Hence, the two ends of the spiral will never coincide and 
the resultant field will never have zero magnitude. 

We have now proved regardless of the physical open 
fraction, there are never missing orders in atom diffrac- 
tion from a material structure unless C3 = 0. In addition, 
the diffraction envelope is no longer described by a sine 2 
function as would be the case if the diffraction were de- 
scribed by a real-valued effective open fraction. 

Our phasor diagram analysis also provides a method 
to bound the error on numerically computed amplitudes. 
If the limits of integration in Equation [2] are replaced by 
±(w/2 — e), the maximum error in resultant amplitude is 
given by the radius of a circle with the curvature of the 
phasor diagram at £ = (w/2 — e), i.e. error in ^ n is less 
than (ty inc /d)R where R^ 1 = d(j> n /d£, | ?=(u) / 2 _ £ ). 

To confirm this theoretical description, we present 
atom diffraction data using a sodium atom beam and 
a silicon nitride grating with a period of 100 nm. We call 
rotation about a grating bar twist as shown in Figure 
Normal incidence defines zero twist. Diffraction pat- 
terns shown in Figure 0] were obtained with the grating 
held at a twist of 0, 12, and then 22 degrees, by scanning 
the position of the hot wire detector. Twist foreshortens 
the grating period and therefore slightly increases the 
diffraction angle. However, the +n and —n diffraction 
orders are nearly equally deflected because the atomic 
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FIG. 3: Experimental geometry. 

de Broglie wavelength is small compared to the grating 
spatial period. 

The feature explored in this study is the relative in- 
tensity in each order, which changes considerably as the 
grating is twisted. Atom flux diffracted into the zeroth 
order decreases when the grating is twisted, but flux 
into the second and third orders first increases. Fig- 
ure 01 shows that the second order intensity is largest 
for an intermediate twist. Variation in the relative in- 
tensity among diffraction orders is expected because the 
projection of the grating viewed from the incident atom 
beam changes with twist. However, a model based on ab- 
sorptive atom optics is not sufficient to explain our data. 
Phase shifts due to van der Waals interactions must be 
included, as discussed earlier. 
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FIG. 4: Diffraction scans with different grating twist. Diffrac- 
tion of both sodium atoms and sodium dimers is visible. First 
order atom diffraction, Ii is located at ±0.3 mm from the th 
order. 

We also present data obtained while scanning the grat- 
ing twist and leaving the detector at one position (Figure 
|3J). In the latter experiment we have measured /„ = \^n\ 2 
for each n=[0,4] while continuously changing the pro- 
jected open fraction. This technique is well suited to 
the task of searching for missing orders, because the pro- 
jected open fraction can be scanned through i, i, j, i.e. 
values that would make missing orders according to the 
optical theory. 

For a perfectly thin grating, twist would not affect the 
open fraction. For our gratings, with a geometry shown 
in Figure El twist does in fact modify the projected open 
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FIG. 5: Data (circles) and models (dashed and solid lines) 
of the intensity in each diffraction order as a function of grat- 
ing twist. The model parameters are: d=I00nm, w=67nm, 
i=II6nm, and a = 3.5°. For dashed lines Cs=0 (using equa- 
tion [^Q, an d f° r solid lines C3—5 eVA 3 (using equation 113^ . 
Statistical error bars for each data point are smaller than the 
circles. 

fraction. Furthermore, due to trapezoidal grating bars, 
van der Waals phase shifts must be carefully analyzed. 
With reference to Figure twisting the grating by angle 
(3 causes slots to appear narrower so the resulting open 
fraction is: 

m = \ ;a<\(3\<p max (10) 

where w is the slot width viewed at normal incidence, 
I is the thickness of the grating, a is the wedge angle of 
the bars, (3 is the twist, (3 max is the maximum twist at 
which any flux is transmitted, and d is the grating period 
viewed at normal incidence. The intensity in different 
diffraction orders then depends on twist as: 

I n (P) =r ] (l3fsmc 2 {ri{P)n) (If) 
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This model without an allowance for atom-surface in- 
teractions was used to predict intensity in the first five 
orders as a function of twist shown in Figure [S] (dashed 
lines), and compares poorly with the data. 




FIG. 6: Top view of the atom beam passing through the 
grating slots. The grating twist is denoted by (3, the wedge 
angle of each bar is a, the grating period is d, the thickness 
is I, and physical open width is w. 

Two features of this model are familiar from standard 
diffraction theory. First, in the limit of small open frac- 
tion the intensity of each order becomes equal. This hap- 
pens at large twist. Second, the orders are completely 
suppressed at certain angles for which n = m7y(/3) , i.e. 
missing orders are predicted. However, the atom beam 
flux is never entirely suppressed until the projected open 
fraction is zero, i.e. there are no missing orders in the 
data. 

Several of the grating geometry parameters are known 
from the manufacturing process and SEM images 0. 
The period is d — 100 nm and the window size is w = 
67±5 nm, so at normal incidence the open fraction is ap- 
proximately 67%. The only physical grating parameters 
left in the optical model are the thickness, I, and wedge 
angle, a. These are constrained by measuring the max- 
imum twist at which any flux is transmitted, and must 
satisfy the condition w = l(ta,n[3 rnax — tana). The data 
show Pmax = 31 ± 2°. With this constraint, there is only 
one free parameter in the optical model, i.e. a. 

The relative intensity for each order is determined in 
the optical model by Ea llll and can be compared to the 
data that were recorded without changing the detector 
gain or incident beam flux. However, the velocity distri- 
bution of the beam broadens the higher order diffraction 
peaks. Each diffraction order has a HWHM in the detec- 
tor plane given approximately by a\ = (Tq + (QnL 2 ^) 2 
where cto is the HWHM of the zeroth order, L is the 
distance between the grating and the detector, ^ is the 
ratio of the HWHM spread in velocity to the average ve- 
locity, and 0„ is the diffraction angle 2]. The velocity 
ratio is 1/15 and the average velocity is 1000 m/s. To al- 



low for the velocity distribution, the relative intensity of 
each diffraction order is multiplied by a~ in the model. 
We have not accounted for the change in 6„ with twist, 
that has the effect of further reducing the recorded in- 
tensity of the higher orders at large twist. 

To include the van der Waals interactions with the 
twisted grating, the model must take account for the 
varying distance to the surface as atoms pass through 
the slots as shown in Figure El i.e. the transverse co- 
ordinate in the grating now depends on the longitudinal 
position £ — > £o + £( z )- The potential due to each interior 
wall is approximated by the van der Waals potential for 
an infinite surface when the atoms are inside the grating 
slots, and zero elsewhere. Then the phase shift due to 
one wall of a twisted grating is: 
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where 9 = [3 ± a with the sign depending on whether 
the classical paths get closer or farther from the grating 
wall as a function of z. Changing the sign of w describes 
the phase shift from the other wall. As before, the field 
amplitude in n th order diffraction is given by an integral, 
but now the limits of integration depend on grating twist 
as shown in Eauation llOl as does <f> v dw given by Equation 

El 

tf n = / e H+n(e)+*v**rW)d£. (13) 

Equation El can now be used to describe intensity in 
each order as a function of van der Waals coefficient, atom 
velocity, and grating twist: I n (Cs,v, /3) = l^n] 2 - When 
C3 is zero, the expression for the intensities reduces to 
Equation ^2 ln comparison, when C3 is not zero the 
model predicts no missing orders and agrees qualitatively 
with the data in Figure (this model is shown with solid 
lines). We note the van der Waals interaction diverts flux 
from the zeroth order into the higher orders, and tends 
to smooth the features given by Equation ^2 Even the 
slight asymmetry in intensity as a function of twist is 
reproduced. To our knowledge this is the first hint of a 
fabricated structure acting as a blazed grating for atom 
waves. 

When used to measure the strength of C3 for sodium 
atoms and a grating made of silicon nitride, we determine 
a value for C3 = St^eVA 3 . Further work on the precise 
shape of the grating, and the van der Waals potential in 
all space due to the structure is needed to reduce this 
uncertainty. 

The maximum asymmetry we observe in the first or- 
der occurs at a twist of ±5° and is I^ 5 =1.1. In 
simulations with a larger wedge angle and larger C3, the 
asymmetry can be as large as 1.5. An asymmetric dis- 
tribution of intensity between the +l st and —1 st orders, 



5 



as this implies, would be useful for atom interferometers 
that only employ the th and +l st diffraction orders. 

Even at normal incidence the van der Waals interac- 
tion has reduced the intensity of the zeroth order diffrac- 
tion by a factor of 0.65 compared to the zeroth order 
flux predicted with C3 = 0. As an extrapolation, if we 
could obtain similar gratings with a 20 nm period, the 
flux transmitted into the zeroth order will be reduced to 
0.06 of the flux predicted with C3 = 0. Hence, van der 
Waals interactions pose a significant problem for sub-100- 
nanometer scale atom optics. 

In conclusion, a novel way to measure the atom-surface 
interaction potential was presented. By twisting a 100 
nm period diffraction grating, we show that atom-surface 
interactions prevent missing orders, and cause asymmet- 
ric diffraction patterns. Both observations are explained 
by a complex transmission function and a phasor analysis 
similar to the Cornu spiral. 

We are indebted to T. Savas and H.I. Smith for fabri- 
cation of the 100 nanometer period material gratings 0- 
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